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ABSTRACT 

We present an updated and improved M^^-a diagram containing 64 galaxies for 
which Mbh measurements (not just upper hmits) are available. Due to new and in- 
creased black hole masses at the high-mass end, and a better representation of barred 
galaxies at the low-mass end, the "classical" (all morphological type) Mbh-c rela- 
tion for predicting black hole masses is \og{Mhh/MQ) = (8.13 ± 0.05) -I- (5.13 ± 
0.34) log[a/200kms~"'^], with an r.m.s. scatter of 0.43 dex. Modifying the regres- 
sion analysis to correct for a hitherto over-looked sample bias in which black holes 
with masses < IO^Mq are not (yet) detectable, the relation steepens further to give 
log(Mbh/M0) = (8.15±0.06)-h(5.95±0.44)log[o-/200kms-^]. We have also updated 
the "barless" and "elliptical-only" Mhh-u relations introduced by Graham and Hu in 
2008 due to the offset nature of barred galaxies. These relations have a total scatter 
as low as 0.34 dex and currently define the upper envelope of points in the Mbh-c 
diagram. They also have a slope consistent with a value 5, in agreement with the 
prediction by Silk & Rees based on feedback from massive black holes in bulges built 
by monolithic-collapse. 

Using updated virial products and velocity dispersions from 28 active galactic 
nuclei, we determine that the optimal scaling factor / — which brings their virial 



products in line with the 64 directly measured black hole masses — is 2. 



J -(-0.7 
-0.5- 



This 



is roughly half the value reported by Onken et al. and Woo ct al., and consequently 
halves the mass estimates of most high-redshift quasars. Given that barred galaxies 
are, on average, located ~0.5 dex below the "barless" and "elliptical-only" Mbh-c 
relations, we have explored the results after separating the samples into barred and 
non-barred galaxies, and we have also developed a preliminary corrective term to the 
velocity dispersion based on bar dynamics. In addition, given the recently recognised 
coexistence of massive black holes and nuclear star clusters, we present the first ever 
(Mbh + Mnc)-(^ diagram and begin to explore how galaxies shift from their former 
location in the A/bh-c diagram. 

Key words: Astronomical Data bases: catalogues — black hole physics — galaxies: 
active — galaxies: nuclei — galaxies: Seyfert — galaxies: quasars: general 



1 INTRODUCTION 

The Mbh-o- relation (Ferrarese & Merritt 2000; Gebhardt et 
al. 2000a) is important because: (i) it enables one to predict 
supermassive black hole (SMBH) masses, Mbh, in galaxies 
for which only the bulge velocity dispersion, a, is known; 



(ii) it allows one to calibrate other relations which can then 
be used to predict SMBH masses in active galactic nuclei 
(AGN) for which a can not be readily measured; and (iii) 
it points toward a physical connection between the nuclei 
of galaxies and the properties of their host bulge. SMBH 
masses themselves, plus their demographics, accretion and 
activity is important for a number of reasons, in particular 
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the influence that SMBHs are thought to have in dictating 
the growth of galaxies and galaxy clusters. 

Most AGN are too distant to spatially resolve any ma- 
terial which is predominantly under the dynamical influence 
of their black hole. Therefore, less direct methods to deter- 
mine the masses of these black holes are required. Rever- 
beration mapping (RM: e.g. Bahcall, Kozlovsky & Salpeter 
1972; Blandford & McKee 1982; Netzer & Peterson 1997) is 
the name given to observations which measure the time de- 
lay between direct continuum emission from a central AGN 
and the echoed emission-line signal from gas clouds in the 
so-called broad line region (BLR: Seyfert 1943) surround- 
ing the AGN (Shields 1974, see also Souffrin 1968). Given 
the constant speed of light, such time delays correspond to 
a distance, r. Coupled with the Doppler-broadened width 
AV^ of the emission lines from the clouds, and assuming 
that their motion is virialised and dominated by the cen- 
tral black hole's gravity (e.g. Gaskell 1988, 2009a; Koratkar 
& Gaskell 1991; Wandel, Peterson & Malkan 1999; Onken 
& Peterson 2002), one can compute the virial product VP 
= rAF'^/G. To convert these VPs into black hole masses 
requires the multiplication by a scaling factor / which is 
related to the geometry and orientation of the clouds and 
effectively converts the measured velocity widths into an in- 
trinsic Keplerian velocity (Peterson & Wandel 2000; Onken 
et al. 2004). 

In diagrams that plot directly measured, and thus hope- 
fully reliable, SMBH masses Mbh obtained from nearby, pre- 
dominantly inactive galaxies versus some other galaxy prop- 
erty such as velocity dispersion, an empirical calibration of 
the above /-factor can be performed by finding the value of / 
which yields the optimal overlapping agreement between the 
virial products and the directly measured black hole masses. 
Adopting a fixed value of / = 3 (Netzer 1990), this cal- 
ibration of RM masses was first explored by Gebhardt et 
al. (2000b) and Ferrarese et al. (2001) using stellar velocity 
dispersions for seven and six AGNs, respectively^ Figure [T] 
presents the results from Onken et al.'s (2004) first ever 
empirical calibration of /, which assumed that the 16 local 
AGNs they studied lie on the inactive galaxy Mbh — cr re- 
lation. The initial Mbh-o" relations for inactive galaxies pos- 
sessed a small total root mean square (r.m.s.) scatter in the 
log Mhh direction of ^^0.34 dex (Merritt & Ferrarese 2001a; 
Tremaine et al. 2002) and therefore appeared well-suited for 
the above task of calibrating the /-factor. However, as more 
(barred) galaxies have been added to the Mbh-o" diagram, 
the r.m.s. scatter has increased above 0.4 dex, particularly 
at the low mass end (e.g. Graham 2008b; Gaskell 2009c). 
Indeed, Graham (2007a, 2008a), Hu (2008) and Graham & 
Li (2009) have revealed that excluding barred galaxies, or 
using only elliptical galaxies, results in the recovery of a 
tight barless, or elliptical-only, Mbh-o" relation, a result re- 
iterated by Giiltekin et al. (2009b). The reason for this is 



^ The frequently used value of 3 is derived by taking the isotropic 
assumption of (T3J3 = voctxd . However, because Gebhardt et al. 
(2000b) and Ferrarese et al. (2001) made use of the emission line's 
full width at half maximum (FWHM), which Netzer took as 2criD , 
they quote an /-factor of 3/4 rather than 3. The RM results which 
we discuss adopt the second moment of the line profile ((Tiino) 
as the measure of the emission line width, and the /-factors we 
describe arc appropriate to such data. 
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Figure 1. Open squares represent 31 predominantly inactive 
galaxies with direct SMBH mass measurements as tabulated by 
Tremaine et al. (2002). The solid (and dashed) line is the rela- 
tion (and associated la uncertainty) reported by Tremaine et al. 
(2002). The shaded region expands this domain vertically by 0.34 
dex — the r.m.s. scatter about the relation. The 16 filled squares 
represent AGN with (reverberation mapping)-derived virial prod- 
ucts from Onken et al. (2004, their Table 3). While in this figure 
an /-factor of 1 has been used to plot these latter points, Onken 
et al. determined that an /-factor of 5.5± 1.7 was required to con- 
vert/elevate these AGN's virial products into black hole masses 
that agreed with the inactive galaxy sample in this diagram (see 
also Woo et al. 2010 who report / = 5.2 ± 1.2). 



not yet clear, but may be partly due to elongated orbits in 
bars or because bars can puff up the (face-on) central veloc- 
ity dispersion as they evolve (see Gadotti & de Souza 2005; 
Gadotti & Kauffmann 2009; Perez et al. 2009; Saha, Tseng 
& Taam 2010). 

Although bars are thin when they first form, vertical 
resonances (Combes & Sanders 1981; Combes et al. 1990; 
Raha et al. 1991) and possibly also hose instabilities (Merritt 
& Sellwood 1994) result in their growth out of the disc plane 
leading to boxy-peanut-shaped bulges (Combes & Sanders 
1981; see also lUingworth 1981 and Kormendy & Illingworth 
1982). Similarly, the torques due to the bar push gas in- 
wards and can result in the formation of pseudo-bulges (Ko- 
rmendy & Kennicutt 2004, Athanassoula 2005). Athanas- 
soula & Misiriotis (2002, their figure 13) have shown how 
the velocity dispersion can increase dramatically, both out 
of the disc plane and along the length of the bar. Therefore, 
barred galaxies have a mechanism by which they can mi- 
grate off the Mbh-cr relation defined by non-barred galaxies. 
Failing to account for this offset population can not only bias 
one's estimates of SMBH masses in AGN, but also impact 
on various evolutionary studies which may be using galaxy 
samples with varying barred galaxy fractions at different 
redshifts. 

In Section 2.1 we introduce the (predominantly quies- 
cent) galaxy sample for which 64 direct supermassive black 
hole mass measurements and velocity dispersions are avail- 
able. These galaxies are used to define updated Mbh-cr re- 
lations. In Section 2.2 we introduce a sample of 28 AGN 
with available velocity dispersions and explain how their 
RM measurements and line widths have been converted 
into virial products. In Section 3 we derive the optimal /- 
factor which brings the AGN in line with the distribution 
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of (MbhjCr) data for the galaxies with direct SMBH mass 
measurements. This factor provides the (sample average) 
calibration needed to convert virial products into SMBH 
masses. We also identify a potential sample selection bias 
which, if uncorrected, results in underestimation of the slope 
of the Mbh-o" relation. A lengthy discussion is provided in 
Section 4, addressing the slope and scatter of the Mbh-o" re- 
lation, sources of uncertainty on measures of Mbh and o, 
and implications of the new /-factor for AGN masses. Fur- 
thermore, a new equation providing a first order correction 
for the influence of bar dynamics on the host galaxy velocity 
dispersion is presented. We also probe how the coexistence of 
SMBHs and nuclear star clusters may alter our understand- 
ing of the Afbh-o" diagram/relation. Our main conclusions 
are summarised in Section 5. 



2 AN UPDATED GALAXY DATA SET, AND 
THE A/bh-ct diagram/relation 

2.1 Inactive galaxies 

Graham (2008b) provided a catalogue of (Mbh,o") values 
for 76 predominantly inactive galaxies having direct SMBH 
mass measurements which were carefully adjusted according 
to their adopted distance. Table 1 from that paper (based 
in part on Hu's 2008 compilation and Ferrarese & Ford 
2005) provides "reliable" entries for 50 galaxies, 36 of which 
are non-barred systems. Giiltekin et al. (2009b) increased 
this samplqfl with the addition of 4 galaxies (Giiltekin et 
al. 2009ajj and we include these here after adjusting their 
masses to the (Hubble const ant) -independent distances re- 
ported by Tonry et al. (2001). Here we have, however up- 
dated the distance moduli that were reported by Tonry et 
al. (2001, their Table 1), and used in Graham (2008b), by 
decreasing their values by 0.06 mag, thereby reducing the as- 
sociated galaxy distances by ~3 per cent, and thus reducing 
the SMBH masses by this same amount. This small correc- 
tion stems from Blakeslee et al.'s (2002, their Section 4.6) 
recalibration of the surface brightness fluctuation method 
using the final Cepheid distances given by Freedman et al. 
(2001, with the metallicity correction). 

In addition to the following 8 galaxies (IC 2560, NGC: 
2974; 3079; 3414; 4552; 4621; 5813; 5846) that were not in- 
cluded by Giiltekin et al. (2009b), we further expand our 
(Mbh,o") catalogue from Graham (2008b, his Table 1) with 
the inclusion of another ten galaxies (see Table [ij. The 
following four were previously considered to have uncer- 
tain SMBH masses (NGC 1068, Lodato & Berlin 2003; 
NGC 3393, Kondratko et al. 2008; Abell 1836-BCG and 
Abell 3565-BCG, Dalla Bonta et al. 2009J3 while the fol- 
lowing six are new galaxies: NGC 253 (Rodn'guez-Rico et 
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Figure 2. Updated Mbh-o" relation, containing 64 galaxies with 
reliable SMBH masses (see Table 1 from Graham 2008b, and Ta- 
ble 1 from this paper). The solid line is derived from a regression 
of log Mbh on logfT, assuming a 10 per cent uncertainty on the 
velocity dispersions a (see equation [TJ • The dashed lines trace 
the 1-sigma uncertainty on this relation, while the shaded area 
extends this boundary by 0.43 dex (the r.ni.s. scatter about the 
relation) in the log Mbh direction. The barred galaxies can be 
seen to dominate the distribution of points below the best-fitting 
line at low-masses. 



al. 2006); NGC 524 and NGC 2549 (Krajnovic et al. 2009); 
NGC 1316 (Fornax A: Nowak et al. 2008); NGC 3368 and 
NGC 3489: (Nowak et al. 2010); This gives a total sample 
of 64 galaxies with reliable SMBH masses. 

Aside from this expansion from 50 to 64 galaxies, 
we have updated the mass for the Milky Way's SMBH 
(Gillessen et al. 2009) and roughly doubled the SMBH 
masses of NGC 3379, NGC 4486 and NGC 4649 (see Ta- 
ble[l} . Recently published masses are increasingly secure due 
to refinements such as the use of triaxial orbit-based models 
rather than spherical or axisymmetric models, and better 
accounting for the range of orbital anisotropics and the in- 
fluence of dark matter (Gebhardt & Thomas 2009; Shen & 
Gebhardt 2009; van den Bosch & de Zeeuw 2010). Follow- 
ing Giiltekin et al.'s (2009b) identification of an error in the 
SMBH masses reported by Gebhardt et al. (2003), we have 
also increased the black hole masses, and their associated 
uncertainties, by a factor of 1.099 in the following 9 galax- 
ies: NGC 821; NGC 2778; NGC 3384; NGC 3608; NGC 4291; 
NGC 4473; NGC 4564; NGC 4697; and NGC 5845. The ex- 
panded Afbh-c" diagram can bee seen in Figure (2] 



^ It is worth noting that there are differences among the adopted 
distances, and thus the adopted SMBH masses, between Graham 
(2008b) and Giiltekin et al. (2009b), an issue already addressed 
by Graham & Driver (2007b, their Section 2.1). 
^ Giiltekin et al. (2009b) also includes 17 galaxies, not tabulated 
in Graham (2008b), for which only upper limits to their SMBH 
masses are available. An additional 105 upper limits can be found 
in Beifiori et al. (2009). 

* The references provided explain why these points arc now con- 
sidered reliable. 



2.1.1 The Mhh-cy relaUon(s) 

For predicting SMBH masses in other galaxies, one obviously 
desires a relation with the minimal amount of scatter in the 
vertical (logAfbh) direction. This is achieved with a non- 
symmetrical ordinary least squares regression of log A'fbh on 
logo- (Feigelson & Babu 1992). Using the BCES code from 
Akritas & Bershady (1996) and assigning a 10 per cent un- 
certainty to the velocity dispersions of the 64 galaxies with 
direct SMBH mass measurements, one obtains 
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Table 1. Expansion and update to Table 1 from Graham (2008b). 



Gal. Id. 



1 



Type 



Dist. 
Mpc 
3 



km s 

4 



Mbh 
10* M0 
5 



New inclusions 



Abell 
Abell 
NGC 
NGC 
NGC 
NGC 
NGC 
NGC 
NGC 
NGC 
NGC 
NGC 
NGC 
NGC 



1836 

3565 

253 

524 

1068 

1316 

2549 

3368 

3393 

3489 

3585 

3607 

4026 

5576 



Milky Way 
NGC 821 
NGC 2778 
NGC 3379 
NGC 3384 
NGC 3608 
NGC 4291 
NGC 4473 
NGC 4486 
NGC 4564 
NGC 4649 
NGC 4697 
NGC 5128 
NGC 5845 



BCG 

BCG 

SBc 

SO 

Sb 

SBO 

SBO [12] 

SBab 

SBab 

SBO 

SO 

SO 

SO 

E3 

SBbc 

E 

SBO 

E 

SBO 

E2 

E2 

E5 

EO 

SO 

El 

E4 

SO 

E3 



157 [la] 

40.7 [la] 

3.5 [2] 

23.3 

15.2 [lb] 

18.6 [3] 

12.3 

10.1 

55.2 [lb] 

11.7 

19.5 

22.2 

13.2 

24.8 



309 [6] 

335 [7] 

109 [8] 

253 

165 [9] 

226 

144 

128 

197 

105 

206 

224 

178 

171 



39^ 



h4 
-5 
11+2 

1+0.1 

-O.OE 
,+2.7 

-1.3 
fO.003 



[6] 
[6] 

O-Haos [11 
8.3+^'I fl2' 
0.084 



1.50 



0.003 [l-^l 
0.80 [14] 



+0.75 



.1+0.02 



0.14^-- [12] 
0.073l°-°i^ [15] 

'^-To2 [16] 
0.058l«:S [15] 

-1.4 



0.34"f 



3.1 
1.3 
1 

1 



Updated data 
0.008 [4] 100 [10] 



0.6 
+0.5 

0.5 
+0.6 
-0.3 
+0.3 
-0.4 



17] 
[17] 
[17] 

[17] 



23.4 

22.3 

10.3 

11.3 

22.3 

25.5 

15.3 

15.6 

14.6 

16.4 

11.4 

3.8 [5] 

25.2 



200 
162 
209 
148 
192 
285 
179 
334 
157 
335 
171 
120 
238 



0.043tHr4 [4] 
0.39+0-26 [18] 

O.lStg-;^ [18] 

4.0l}-0 [19] 

0.17+0;°^ [18] 

2.0lS:^ [18] 

3.3l«:^ [18] 

1.2+0;^ [18] 

56l* [20] 

O.eotHg [18] 
4711° [21] 
1.8l°:? [18] 
0.45t°-}^ [22 
2.6+?-f - - 

"1.5 



[18] 



Unless otherwise specified, the distances have come from Tonry 
et al. (2001), after reducing their distance moduli by 0.06 mag 
(see Section |2. 11 1. This small adjustment has been applied to all 
the galaxies from Graham (2008b, his Table 1) which used the 
Tonry et al. (2001) distance moduli. Unless otherwise specified, 
the velocity dispersions, a, are the weighted values from Hyper- 
LedE0 (Paturel et al. 2003) as of March 2010. The SMBH masses, 
Mbh, have been adjusted to the distances given in column 3. 
References: la = NED; (Virgo -|- GA -|- Shapley)-corrected Hub- 
ble flow distance of the BCG's host cluster; lb = NED: (Virgo -|- 
GA -|- Shapley)-corrected Hubble flow distance; 2 = Rckola et al. 
(2005); 3 = Madore et al. (1999); 4 = Gillessen et al. (2009); 5 
= Karachentsev et al. (2007); 6 = Dalla Bonta et al. (2009); 7 = 
Smith et al. (2000); 8 = Oliva et al. (1995); 9 = Nelson & Whittle 
(1995); 10 = Merritt & Ferrarese (2001a); 11 = Rodriguez-Rico 
et al. (2006), a factor of 2 uncertainty has been assigned here; 12 
= Krajnovic et al. (2009); 13 = Lodato & Bertin (2003); 14 = 
Nowak et al. (2008); 15 = Nowak et al. (2010); 16 = Kondratko 
et al. (2008); 17 = Giiltekin et al. (2009a); 18 = Gebhardt ct al. 
(2003), Giiltekin et al. (2009b); 19 = van den Bosch & de Zeeuw 
(2010); 20 = Gebhardt & Thomas (2009); 21 = Shen & Gebhardt 
(2009); 22 = Neumayer (2010). 

log(Mbh/AfQ) = (8.13±0.05)+(5.13±0.34) log[cr/200kms"^], (1) 

with a total r.m.s. scatter of A = 0.43 dex, and an intrin- 
sic scatter of e = 0.32too4 dex, in the logMbh direction. 
This relation is shown in Figure [2l For reference, the equiv- 
alently produced relation using all 50 galaxies from Gra- 
ham (2008b) has a slope and intercept of 4.87 ± 0.36 and 



8.12 ± 0.06, respectively, while Hu (2008) reported values 
of 4.59 ± 0.32 and 8.14 ± 0.06. Using our updated data set, 
the slope drops to 4.86 ± 0.31, while the intercept remains 
unchanged, when an uncertainty of 5 per cent — the value 
used by Tremaine et al. (2002) and Giiltekin et al. (2009b) 
— is assigned to the velocity dispersions. Using Tremaine 
et al.'s (2002) modified FITEXY routine (Press et al. 1992) 
yields consistent results, with a slope and intercept of 4.86 
and 8.15, respectively (when assigning a 5 per cent uncer- 
tainty to the velocity dispersions). Of the two commonly 
used uncertanties within the literature (namely, 5 and 10 
per cent), we have elected to proceed using a 10 per cent 
uncertainty on our velocity dispersions, as an inspection of 
HyperLeda (Paturel et al. 2003, see also Section 6 in Nowak 
et al. 2010) reveals that a 5 per cent uncertainty is proba- 
bly overly optimistic. Ideally though, one would like to have 
more accurate measurements and measurement errors for 
the velocity dispersions (see Section [4. 2p . 

When using (only) the 44 non-barred galaxies, and a 10 
per cent uncertainty on the velocity dispersions, the above 
regression of log Afbh on log a gives the relation 

log(Mbh/M0) = (8.25±0.06)+(4.57±0.35) log[cr/200kms"^], (2) 

with a total r.m.s. scatter of A = 0.37 dex, and an intrinsic 
scatter of e = Q.2<d't°o% <is^' i^i ^^^ logMbh direction0 As 
already noted by Graham (2008a) and Hu (2008), a similar 
relation is obtained when using only the elliptical galaxies 
(see Table[5]). The r.m.s. scatter of the elliptical-only Mbh-c 
relation is 0.34 dex, notably less than the value of 0.43 dex 
for the standard (full sample) Mbh-o" relation. 

When using (only) the 20 barred galaxies, and a 10 per 
cent uncertainty on the velocity dispersions, the regression 
of log Mbh on log a gives 

log(Mbh/M0) = (7.80±0.10)+(4.34±0.56) log[cr/200kms"^], (3) 

with a total r.m.s. scatter of A = 0.36 dex. For convenience, 
all of these Mbh-o" relations are tabulated in Table [^j This 
barred Mbh-o" relation is 0.45 dex (at a = 200 km s~^) below 
the relation defined by non-barred galaxies, and reiterates 
the offset nature of the two populations first noted by Gra- 
ham (2007a, 2008a,b) and Hu (2008). 

Whilst we adopt the above approach for this section, 
we do note that Section 13.11 points out, for the first time, 
a potential sample bias which, once corrected for, results 
in steeper slopes for the various Mbh-o" relations. These are 
given in the second half of Table[2l In Section Fl.ll we discuss 
the causes for the change in slope of the Mbh-o" relation. 



2.1.2 Predicting M^^ 

When predicting the SMBH mass of a new galaxy for which 
one knows the velocity dispersion, the associated maximum 
1-sigma uncertainty on the black hole mass — acquired by 
assuming uncorrelated errors on the velocity dispersion a 
and both the slope and intercept of the Mbh-o" relation — 
can be determined using Gaussian error propagation. For 



° When an uncertainty of 5 per cent is assigned to the velocity 
dispersion, the slope of the "barless" Mbh-o" relation drops to 
4.32 it 0.34, while the intercept basically remains the same. 
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Table 2. Assorted log{Mbh/M0) 
relations. 



a + /31og((T/200 km s"!) 



CT- X crfinJG, 



(5) 



Sample (size) 


a 


13 A 


log A/bh 
[dcx] 


^intrinsic 

[dcx] 


BCES 


regression of 


log A/bh on log o- 






Full (64) 


8.13 ±0.05 


5.13 ±0.34 


0.43 


n nO+0.06 


Barred (20) 


7.80 ±0.10 


4.34 ± 0.56 


0.36 


n 97+0.09 


Non-barred (44) 


8.25 ±0.06 


4.57 ±0.35 


0.37 


29+°°'^ 


Elliptical (25) 


8.27 ±0.06 


4.43 ±0.57 


0.34 


27+°-"^ 


BCES 


regression of 


log cr on log A/bh 






Full (64) 


8.15 ±0.06 


5.95 ± 0.44 


0.46 


'-'■■'^-0.05 


Barred (20) 


7.94 ±0.18 


5.40 ± 1.22 


0.41 


U-3t_Q ;^Q 


Non-barred (44) 


8.24 ±0.06 


5.32 ± 0.49 


0.41 


S0+°°^ 


Elliptical (25) 


8.22 ±0.09 


5.30 ±0.77 


0.37 


29+°"* 



for each AGN. Modulo the yet-to-be-determined /-factor, 
this virial product represents the mass of each AGN's central 
black hole. 

When multiplying uncorrelated numbers with errors, 
the relative error on the product is the square root of the 
sum of the squares of the relative errors in the individual 
numbers. That is 

(xi ± Sxi) X (xj ± Sxj) X ...{xm ± Sxm) 



— {Ili^l^NXi 



1± 



fSxi 



i:(f 



Based on a 10 per cent uncertainty on the velocity dispersions cr. 
The intrinsic dispersion ^intrinsic pertains to the log A/bh direc- 
tion. 

the linear equation y = {b± 5b){x ± Sx) + (a ± Sa), one has 
an error on y equal to 

Sy = ^J{dy/dbY{ShY + {dy/da)^{Say + {dy/dx)'2{Sx)^ 

= y'x^iSb)^ + (Sa)^ + b'2i6xy. 

In the presence of intrinsic scatter in the j/-direction, denoted 
by e, the uncertainty on y is 



Sy^ v/x2(5fe)2 + (5a)2 + fe2(5a:)2 + e2. 

Using the standard Mbh-o" relation, given by equation [T] we 
have that x = log(o'/200kms~^), so dx/da — l/[ln(10)(7], 
and therefore 

{SlogMbb/MQf = [log(cr/200 km s-')]^ (0.34)^(0.05)^ 
+ [5.13/ ln(10)]^[<5a/o-]^ + (0.32)^ (4) 

This represents the uncertainty when predicting a new black 
hole mass from a velocity dispersion measurement a iL 5a. 
Of course, if one knows the morphological type of the galaxy 
in question, then the relevant equation from Table [2] can be 
used to predict a more accurate black hole mass. 

Before proceeding, we note that several ways of dis- 
tinguishing pseudo bulges from classical bulges have been 
proposed in the literature. This includes morphology, kine- 
matic properties, the Sersic index and the distance from the 
Kormendy (1977) relation. Unfortunately there is no good 
agreement between these methods and often authors using 
one method criticize the method(s) used by others. On the 
other hand, whether a galaxy has a bar, or not, is a more 
clear cut question. We thus believe that this is a much safer 
way of addressing departures from the Mbh-o" relation. 



2.2 Active galcixies 

We have started with the homogenized time delays, r, and 
broad line region velocity dispersion measurements for the 
35 AGN listed in Table 6 of Peterson et al. (2004). Errors on 
both measurements are available, as are multiple measure- 
ments for many galaxies[j Below we describe how this data 
has been combined to acquire the single virial product. 



^ Potential problems with single epoch data are described in 
Gaskell (2009b). 



.When the numbers are correlated, the relative errors are 
simply added. Therefore, in deriving the virial product, vp — 
CTcontcjinc/G', froiu the Tcont and (Tiinc valucs tabulated by 
Peterson et al. (2004), one has a relative uncertainty given 
by 



Svp 
vp 



( 5(Tlinc , <5(Jlinc 



V Clinc 



0"lino 



(6) 



These values are provided, for multiple measurements, in 
the final column of Table 6 from Peterson et al. (2004) and 
adopted here. 

Assuming that these individual measurements, vpi ± 
5vpi, of the virial products are distributed normally about 
the true value VP, one can use a maximum likelihood 
analysis to determine the (error-weighted average) optimal 
value and its associated uncertaintyU From the treatment 
of weighted averages by Taylor (1997, Chapter 7), one has 
that 



VP = 



2^1=1 



JV ^i'"P^ , 1 

— , where Wi = 



and the uncertainty on the value of VP is given by 
5VP^ ^ 



2^i = l, 



(7) 



(8) 



The resulting (single) virial product for each galaxy is listed 
in Table |l(c.f. Peterson et al. 2004, their Table 8; and Woo 
et al. 2010, which appeared while we were preparing this 
work). Given that we ultimately find a different /-factor to 
Woo et al. (2010), we felt that it was beneficial to show, in 
Table IH the exact data that we have used. 

We have excluded PG 1211+143 and IC 4329A because 
the associated uncertainty on their virial products give val- 
ues consistent with zero (Peterson et al. 2004). We have 
also excluded the virial product for 3C 390.3 because it has 
a double-peaked emission line profile. While NGG 5548 has 
the most extensive reverberation-mapping of all the AGN 
(Bentz et al. 2007), frustratingly, it also has an irregular 
emission line profile. As noted by Zhu & Zhang (2009) and 
Zhu, Zhang & Tang (2009; see also Wandel, Peterson & 
Malkan 1999), the determination of the line width is prob- 
lematic for these last two galaxies. To err on the side of cau- 
tion, we present our analysis both with and without NGG 



° The one simplification is to average the slightly non- 
symmetrical errors associated with the virial product measure- 
ments, providing the values of Svpi used above. 
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5548, finding consistent results. We have however included 
12 new galaxies (see also Woo et al. 2010) for which velocity 
dispersions were not available at the time of Onken et al.'s 
(2004) analysis. Table |4] lists 30 AGN with reverberation- 
mapping measures and velocity dispersions, two of which 
we exclude, and 13 of which belong to barred galaxies (see 
Bentz et al. 2009b, their Table 5). These data are plotted in 
Figure [31 



3 ANALYSIS: THE F- FACTOR 

Figure [3^ presents an updated version of Figure [l] Again, 
an /-factor of 1 has been applied to the AGN virial prod- 
ucts. Immediately apparent is that the AGN virial products 
no longer appear quite so offset from the directly measured 
SMBH masses, instead, the two populations somewhat over- 
lap. The larger discrepancy seen in Figure[l]appears to have 
been a result of sample selection, specifically, due to the past 
under-representation of barred galaxies with direct SMBH 
masses. 

As noted previously, the /-factor is the normalization 
of the AGN virial product which is used to estimate Afbh 
such that 



Mb 



/x 



G 



(9) 



Here we seek the value of / which best matches the AGN in 
the Mbh-o" diagram with the Mbh-o" relation defined by the 
local sample of galaxies with direct SMBH mass measure- 
ments. 

The x^ value that we minimise to determine this opti- 
mal value of / is given by the expression 



[log(Mbh,»/M0) -a- 61og(cr,/200kms-i)]2 



1 iVPi 



2 + [51og(Mbh,,/M0)]2 



.(10) 



AT 

where a and b are the intercept and slope of the Mbh-o" 
relation, Mbh.i comes from equation |9l 51og(Mbh,i/M0) is 
given in equation [J] and 5VPi is the uncertainty derived us- 
ing equation [8] and listed in Table ID Matching all 28 AGN 
to the standard (full sample) Mbh-o" relation given by equa- 
tion [T] gives / — 3.8lQg. This result is shown in Figure |3t- 
Excluding NGC 5548, / = 3.6l°;^. 

As we have effectively seen, both here and in the liter- 
ature, the inclusion of varying numbers of barred galaxies 
in the Mbh-o" diagram alters the best-fitting Afbh-cr relation 
that one obtains. We have therefore attempted to explore 
a derivation of the /-factor using only barred galaxies, and 
then again using only non-barred galaxies. 

A visual inspection of the barred galaxies in Figure |3j3 
and Uja reveals that the virial products from the 13 barred 
AGN overlap with the territory occupied by the directly 
measured black hole masses of barred galaxies. That is, an 
/-factor closer to a value of unity appears apt, and a formal 
treatment produces a value of / = 2.3lQg. 

Figure [3]:; reveals that 4 of the 15 unbarred AGN have 
black hole masses which agree exceedingly well with the re- 
lation defined by the 44 non-barred quiescent galaxies when 
the AGN virial products are multiplied by an /-factor of 
1. At the same time, 11 of the AGN appear to require an 
/-factor significantly greater than a value of ~l-3. Collec- 
tively, all 15 non-barred AGN generate an optimal /-factor. 
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Figure 3. Panel a) An updated version of Figure [T] using the 
galaxies described in Section 2. One can see that the full sample 
of 28 AGN (filled squares) with (reverberation mapping)-derived 
virial products overlap with some of the 64 galaxies with direct 
SMBH mass measurements (open squares). The lines and shaded 
region pertain to equation [T] In panels b) and c), both AGN and 
galaxies with bars are denoted by a star, while those without by 
a circle. Filled and open symbols pertain to our sample of AGN 
and galaxies with direct SMBH mass measurements, respectively. 
In all panels, as in Figure [T] an /-factor of 1 has been used to 
plot the AGN virial products by the filled symbols. 



relative to equation [2l of — T.Q'^W (/ = 6.51^4 when ex- 
cluding NGC 5548), while using the 11 most offset of these 
15 AGN gives an unlikely /-factor in excess of 10. Such a 
high value of / produces SMBH masses for four of the non- 
barred AGN which are inconsistent with the upepr envelope 
of points in the A/bh-f" diagram defined using direct SMBH 
masses. 

While none of these 11 galaxies display evidence of 
a bar, we do note that three are peculiar in appearance 
(Ark 120, Slavcheva-Mihova & Mihov 2010; Arp 151; and 
3C 120), two may be too distant or small to discern a bar 
(PG 2130-f 099 and Mrk 202, respectively), and another two 
(Mrk 279 and Mrk 590) were modelled with an "inner bulge" 
by Bentz et al. (2009c), a feature that may be related to 
bars (Peng et al. 2002). If we have over-looked the presence 
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Figure 4. Similar to Figure [S] except that the optimal /-factor 
(inset in figure) to bring the AGN virial products in line with 
the corresponding M^ij-o" relation has been used. As in Figure [S] 
the filled symbols in panels a), b) and c) correspond to the full, 
barred and non-barred AGN sample, while the open stars and 
circles correspond to the barred and non-barred galaxies with di- 
rect SMBH mass measurements. A linear regression of log M^y^ on 
log (T for the full, barred and non-barred sample of galaxies with 
direct SMBH mass measurements (see Section l2.1|l has been used 
to construct the three different relations (given by equations [T] [3] 
andlJl shown in panel a), b) and c), respectively. 



of bars in (some of) our allegedly non-barred AGN sample, 
then our separation of AGN galaxies would obviously be in 
error and we would thus be comparing the wrong types of 
galaxies in Figure 2]:, and likely over-estimating the value of 
/. Unfortunately, we feel that it may therefore be more ap- 
propriate at this time to prefer the /-value acquired without 
any attempted division into barred and non-barred galaxies 
(i.e. Figure |4^). 

Since commencing this project, several other potential 
mechanisms which may cause offset behavior in the A/bh-f7 
diagram have come to our attention, all of which move galax- 
ies below or rightward of the upper envelope of points in the 
Mbh-o" diagram. First, as we discuss later in Section [4.41 the 
offset nature of the barred galaxies in the Mbh-o" diagram 
may also be, in part, due to the exclusion of the nuclear 



star cluster mass. If so, this would further negate the appro- 
priateness for the separation of galaxies solely on the basis 
of whether or not they contain a bar because (i) the mass 
of the neglected nuclear cluster may be important and (ii) 
non-barred galaxies can also be nucleated. A second issue, 
briefly raised in section 14.51 pertains to radiation pressure 
from the AGN. This may cause the derived virial product 
to erroneously fall below the Mbh-f" relation. If the broad 
line region has a flattened spatial distribution, then inclined 
AGN may appear to have narrower emission-line widths rel- 
ative to those observed from an edge-on orientation, and 
hence their black hole mass may be under-estimateqj. In 
addition, if some AGN are effectively ignited by tidal inter- 
actions or minor mergers driving gas inward (e.g., Dasyra et 
al. 2007), this external trigger may potentially also elevate 
the velocity dispersion of the galaxy. 

Regrettably, the above three issues effectively handicap 
our ability to proceed as we had hoped, and for the present 
time we fall back to the standard approach used to date, 
which is to neglect the morphological type in the analysis of 
the /-factor. Nonetheless, our increased sample size is more 
representative of the galaxy population at large, and we feel 
that we are able to present a more appropriate measurement 
of this calibration factor. 



3.1 Sample selection bias, and a re-derivation of 
the /-factor 

Collectively, we have not yet managed to directly, from spa- 
tially resolved kinematics, measure the masses of SMBHs 
below about 10*' M0 at the centres of galaxies. This does not 
imply that such "intermediate mass black holes" (IMBHs) 
do not exist, indeed, evidence is accumulating that they 
probably do exist (e.g. Pox 52, Thornton et al. 2008; 
NGC 4395, Filippenko & Ho 2003; Greene & Ho 2004,2007; 
Dong et al. 2007; Naik et al. 2010; Seth et al. 2010). We do 
however note that if some fraction of IMBHs form outside of 
galactic cores, or are kicked out by a gravitational radiation 
recoil event (e.g. Komossa & Merritt 2008, and references 
therein), then the relatively long dynamical friction time- 
scale for them to inspiral to the centre may result in some 
fraction of IMBHs still wondering outside of galaxy cores. 
Indeed, evidence for such a non-central IMBH may have al- 
ready been discovered in galaxy ESO 243-49 (Soria et al. 
2010; Wiersema et al. 2010). 

The non-detection of IMBHs, with resolved gravita- 
tional sphere's of influence, at the centres of galaxies sug- 
gests that a (typically neglected) selection bias may exist 
within the A^bh-o" diagram. If our survey selection is such 
that SMBHs with masses less than r^l(f Mq are excluded, 
then it becomes necessary to modify the type of linear re- 
gression which is used if one is to construct a non-biased 
Mbh-o" relation. To avoid this sample selection bias requires 
a regression which minimises the residuals in the log o direc- 
tion. A discussion of this problem and solution can be found 
in Lynden-Bell et al. (1988. their Figure 10). 

We have therefore repeated the previous analysis using 
an Mbh-o" relation constructed from a linear regression of 



® AGN with BLRs significantly more edge-on than average could 
actually end up with over-estimated SMBH masses 
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Figure 6. Similar to Figure[4l except that a regression of log a on 
log Mbii for the galaxies with direct SMBH mass measurements 
has been used to construct the three linear regressions in each 
panel (see Table [2]l. Excluding NGC 5548 gives / = 2.6^ 
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-0.4' 



2.3_Qg and 4.9^^' in panels a), b) and c), respectively. 



log a on log A/bh for (i) the full sample (see Figure [Sjl , (ii) 
the barred galaxies and (iii) the non-barred galaxies. The 
resulting relations, provided in Table [21 are steeper than 
previously obtained^ Consequently, this results in a reduc- 
tion to the optimal f-factors which are shown in Figure [B] 
Using the 28 AGN and 64 predominantly-inactive galaxies, 
the optimal "bias-free" /-factor is 2.S^fj^. This value is a 
factor of 2 less than reported by Onken et al. (2004) and 
Woo et al. (2010), but in good agreement (perhaps for dif- 
ferent reasons) with the value of S.l^j'g from Marconi et 
al. (2008, see section [53J. Excluding the non-barred AGN 
NGC 5548 resuhs in /-factors of / = 2.Qt°ol, 2.3l^;^ and 
4.9^i' in Figures[6K), b) and c), respectively. 



1" The slope of the VP-cr relation for the 28 AGN is 5.09 ± 0.88, 
consistent with the associated slope of the M^^-a relation (5.95 it 
0.44) for the 64 predominantly inactive galaxies. We have there- 
fore not explored changes in / as a function of black hole mass or 
velocity dispersion. 



4 DISCUSSION 

4.1 The slope, and scatter, of the Mbh-cr relations 

The slope of the Mbh-o" relation has received a lot of at- 
tention since Gebhardt et al. (2000a) reported a value of 
3.75 ±0.3 while Ferrarese & Merritt (2000) reported a value 
of 4.8 ± 0.5. Such slopes had been predicted because of 
the way feedback from accretion onto the central massive 
black hole re-directs energy and momentum back into a 
galaxy, establishing a correlation in which Afbh scales with 
o-"' (Fabian 1999) or ct^ (Silk & Rees 1998; Haehnelt, Natara- 
jan & Rees 1998). However, disagreement over the exact 
slope has continued, with Tremaine et al. (2002) claiming 
a value of 4.01 ± 0.32 and Ferrarese & Ford (2005) advo- 
cating 4.86 ± 0.43. An insightful discussion as to why those 
measured slopes varied, depending on the type of regression 
used, is provided by Novak, Faber & Dekel (2006), see also 
Feigelson & Babu (1992) for background understanding. 

Relative to the "classical" Afbh-o" relation based on the 
use of galaxies of every morphological type, Graham (2007a, 
2008a) and Hu (2008) have revealed that there is a tighter 
relation based on the use of elliptical-only galaxies, or non- 
barred or pseudo bulge galaxies. In addition, these new re- 
lations have a different, shallower, slope than the one ob- 
tained using every galaxy. Using a symmetrical regression 
analysis, Giiltekin et al. (2009b) recently reported a slope of 
4.21±0.45 for 41 allegedly non-barred galaxies with available 
SMBH mass measurements, almost identical to the slope of 
4.28 previously reported by Graham (2008b) when using a 
sample of 36 non-barred galaxies and the same 5 per cent un- 
certainty assigned to the velocity dispersions. Using a non- 
symmetrical regression of log Afbh on log cr, Graham (2008a) 
had also shown that the slope of the "elliptical-only" Afbh-cr 
relation was as low as 3.68±0.25. Hu (2008) also constructed 
such an "elliptical-only A/bh-o" relation", finding a slope of 
3.82 ± 0.36 when using a 5 and 10 per cent uncertainty for 
the early- and late-type galaxies, respectively. Using an un- 
certainty of 5 per cent for the velocity dispersions, Giiltekin 
et al. (2009b) subsequently reported a slope of 3.96 for their 
"elliptical-only" Afbh-o" relation, and a slope of 3.86 when 
using their early- type galaxies. 

We have shown, with an updated and expanded data 
set, that the above slopes have increased. The reason is in 
part due to new data and the doubling of some previous 
SMBH masses at the high-mass end. Table [5] reveals that 
a symmetrical treatment of the data — obtained by aver- 
aging the slopes obtained when regressing log Afbh on log cr 
and logo- on logA^bh — yields a slope of (4.57-f5.32=) 4.95 
and (4.43-f 5.30=) 4.87 for the non-barred and elliptical-only 
galaxies, respectively. While this is equivalent to the (sym- 
metrical regression)-derived slope reported by Ferrarese & 
Merritt (2000) and Ferrarese & Ford (2005), one needs to 
keep in mind that their slope pertained to a classical (barred 
plus unbarred) galaxy sample. Curiously, the optimal value 
for the slope of the barless and elliptical-only Afbh-"" re- 
lations has changed from ~4 to ^-^5, consistent with the 
prediction by Silk & Rees (1998) based on feedback from 
SMBHs in bulges built by monolithic-collapse. Given the po- 
tential sample bias which currently excludes SMBH masses 
less than about one million solar masses, the relations in the 
lower half of Table [2] should be preferred; these are consis- 
tent with a slope of 5. The slope for the full galaxy sample 
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Figure 5. Similar to Figure |2] except that a regression of logcr on log^A/^ij for the 64 galaxies with direct SMBH mass measurements 
has been used (see Table |2}. 



(barred plus unbarred galaxies) is steeper still. The reason 
for this is also partly due to the increase of some SMBH 
masses at the high-mass end, and the inclusion of more 
barred galaxies at the low-mass end. Given that the barred 
gala^xies appear to define their own offset relation relative to 
the non-barred galaxies, coupled with the observation that 
they have velocity dispersions which only span the lower- 
half of the range spanned by the non-barred galaxies, their 
increeised numbers results in a steepening of the "classical" 
Mbh-cr relation. Using all 64 galaxies yields a slope between 
5 and 6 (see Table[2J, a result which appears to be subject to 
the relative numbers of differing morphological type (barred 
versus unbarred galaxies) that one includes. Although, this 
latter remark may need to be revoked depending on the role 
that nuclear star clusters play. 

Finally, it is worth noting that the scatter about the 
"classical" Mbh-o" relation remains in excess of 0.4 dex. It 
would thus appear that the velocity dispersion is not the fun- 
damental/sole parameter driving the host galaxy connection 



with the SMBHsEB The Mbh-n relation (Graham & Driver 
2007a) , which uses the major-axis Sersic index n, has a total 
scatter of only 0.31 dex. In addition, (Graham 2007b) has 
shown that the Mhh-L relation (McLure & Dunlop 2002; 
Erwin et al. 2003; Marconi & Hunt 2003) currently has a 
total scatter of only 0.33 dex when (i) based on near-IR 
CCD images rather than optical photographic plates, (ii) 
using R^'" rather than i?^'''* modelling of the bulge light, 
(iii) after the correct identification of disc galaxies and thus 
the correct separation of bulge and disc fight, and (iv) af- 
ter applying internal dust corrections to spiral and lenticular 
galaxies. When such corrections are made, rather than using 

^^ Technically, it is the intrinsic scatter rather than the total scat- 
ter which reveals what physical quantity may control the SMBH 
mass. However, uncertainty in the measurement errors of the ve- 
locity dispersion render the intrinsic scatter a rather unreliable 
quantity. Nonetheless, using a 5 per cent uncertainty, Giiltekin et 
al. (2009b) report an intrinsic scatter of 0.44 it 0.06 dex for their 
full sample A/bh-o" relation. 
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data from de Vaucouleurs et al. (1991), as done by Giiltekin 
et al. (2009b), the uncertainty on the slope and intercept 
of the Mhh-L relation are reduced, and the intrinsic scatter 
drops from 0.38 dex to 0.30 dex. Consequently, one may be 
tempted to conclude that the bulge luminosity and major- 
axis Sersic index appear to be capable of predicting black 
hole masses with more accuracy than the stellar velocity dis- 
persion. However, we note that Graham & Driver (2007a) 
contains only six barred galaxies, while Graham's (2007b) 
re-analysis of the Marconi & Hunt (2003) galaxy sample con- 
tains only four barred galaxies. If, after an accurate mod- 
elling of the stellar light distribution, barred galaxies are not 
found to be offset in say the Mbh--f'buigo diagram, and the 
scatter remains small, it would suggest that the a values 
may be responsible for these galaxies' offset nature in the 
Mbh-cr diagram. If, on the other hand, barred galaxies are 
similarly offset in the Afbh-i diagram as they are in the Afbh- 
(T diagram, it may be suggestive of an underweight SMBH 
mass relative to the relation / upper-envelope defined thus 
far (see Batcheldor 2010). 



4.2 Sources of uncertainty in Afbh and a 

There are many sources of uncertainty in the values of Afbh 
and (7. While typically not spoken about, they are present 
in all the previous studies of the Afbh-o" relation, and also 
this one. While Peterson (2010) have already provided cau- 
tionary remarks on the reliability of reverberation-mapping 
derived SMBH masses, Cappellari et al. (2010) do the same 
for direct measurements of gas and stellar dynamics around 
SMBHs. Attempting to resolve these issues is not only be- 
yond the intended scope of this paper, but would certainly 
make it quite unwieldy. Nonetheless, we felt it was appropri- 
ate, and hopefully helpful, to list a few of our own concerns. 



4-2.1 Cautionary remarks regarding Afbh 

In regard to the SMBH mass, we list three issues. As we have 
seen from Table [T] observers continue to refine/modify their 
SMBH mass measurements as new data becomes available 
and new techniques are implemented (e.g., Valluri, Merritt 
&: Emsellem 2004; Cappellari et al. 2010). Sometimes this 
reveals an under-estimation of past error bars on SMBH 
masses. Point 1) Using triaxial, rather than oblate or spher- 
ical, models can result in factor of two changes and possibly 
more (van den Bosch & de Zeeuw 2010). 2) Failing to al- 
low for a NC with a different stellar mass-to-light ratio from 
the underlying bulge will bias one's mass measurement of a 
SMBH. 3) A proper dynamical treatment of bars in regard 
to the mass measurement of SMBHs, has in general not been 
undertaken, and Krajnovic et al. (2009) have remarked that 
"The bias introduced by modelling a likely barred galaxy 
using a model with a static, axisymmetric potential has so 
far not been well explored." 



4-2.2 Cautionary remarks regarding a 

At the top end of the Afbh-o" relation, one is working with 
massive elliptical galaxies which can have steep velocity dis- 
persion gradients near their centre. We are referring to the 
profile beyond the SMBH's sphere of influence, where the 



velocity dispersion profile is steep due to the high central 
concentration of stars, as traced by the Sersic index (e.g. 
Trujillo, Graham & Caon 2001). Under poorer seeing condi- 
tions, the measured central velocity dispersion will be under- 
estimated as one effectively samples more of the surrounding 
light which has a lower velocity dispersion (cf. Graham et 
al.'s 1998 value of 353 ± 19 km s~^ and D'Onofrio et al.'s 
1995 value of 420 ±27 km s"^ for NGC 1399). Furthermore, 
given the declining nature of velocity dispersion profiles with 
radius, especially in big elliptical galaxies, the use of an av- 
erage luminosity- weighted velocity dispersion over a larger 
aperture should be systematically different from the cen- 
tral velocity dispersion. The use of total, infinite aperture, 
luminosity-weighted velocity dispersion, rather than central 
velocity dispersion, would result in an Afbh-o"totai relation 
which is considerably steeper than the Afbh-o"contrai relation. 
This is readily appreciated by looking at the aperture ve- 
locity dispersion profiles in Graham & CoUess (1997, their 
figure 8) and Simonneau & Prada (2004, their figure 5). 

At the lower-mass end of the Afbh-c" relation, when deal- 
ing with disc galaxies, the velocity dispersion can be over- 
estimated due to strong rotational gradients within the inner 
region used to measure a. The transition in rotational ve- 
locity (Kot), from say -|-220 km s"^ to -220 km s"^ (the 
rotational velocity of our Sun, e.g. Liu & Zhu 2010) across 
the centre of a disc galaxy can happen over a small radial 
range. Given that some fraction of this rotational velocity 
change will occur within one's adopted aperture or fraction 
of the spectrograph's slit-length, the absorption line profile 
from the integrated flux can be highly broadened in non- 
(face on) discs, thereby significantly increasing the measured 
velocity dispersion above the true value (see Epinat et al.'s 
2010 work with emission lines). Needless to say, given that 
the apparently offset galaxies in the Afbh-o" relation tend to 
be disc galaxies (Graham 2008a; Hu 2008; Giiltekin et al. 
2009b), this scenario is rather interesting. In passing we re- 
mark, perhaps for the first time, that this issue may be a 
factor in explaining the high velocity dispersion measure- 
ments reported for some compact galaxies at high redshifts 
(van Dokkum, Kriek & Franx 2009). 

Velocity dispersions can come from spectrograph slits 
placed along the major axis, the minor axis, or the bar axis. 
Apertures and slits can be of varying radius and length 
(and thus sample varying fractions of a bulge's half-light 
radius); together with different galaxy distances, plus vary- 
ing bar lengths, this can all act to modify the velocity dis- 
persion that is actually measured. Detailed measurements, 
from integral field spectrographs, for the 64 galaxies with 
direct SMBH mass measurements would be a welcome con- 
tribution (see Batcheldor et al. 2005). Such uniform data, 
sampling spectral lines that are not particularly sensitive to 
young stars (such as the Ca Triplet lines) , and whose acqui- 
sition is admittedly beyond the scope of this paper, would be 
preferable to the publicly available, hetereogeneous mix of 
velocity dispersions which are, thankfully, homogenised by 
HyperLeda (Paturel et al. 2003) following the precepts of 
McElroy (1995) to produce the galaxy velocity dispersions 
that we and others use. Having a combination of data from 
many authors does, however, reduce the impact that sys- 
tematic biases — which may be present in some individual 
observing programs — can have. 

The problem with young stars, alluded to above, is that 
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they shine very bright, and can bias, and even dominate, 
one's luminosity weighted velocity dispersion (e.g. Wozniak 
et al. 2003). One may therefore have a situation where a cold 
disc has been able to form new stars, and this young nuclear 
disc dominates the central emission but does not represent 
the dynamics of the classical bulge. Similarly, one may have 
a significant contribution from a relatively cold nuclear or in- 
ner bar. This may result in "sigma drops" where the central 
velocity dispersion is actually depressed (e.g. Franx et al. 
1989; D'Onofrio et al. 1995; Graham et al. 1998; Heraudeau 
& Simien 1998) or the value at the ends of the inner bar is 
lowered (de Lorenzo-Caceres et al. 2008). 

Studying the nuclear region of three barred galaxies, 
Emsellem et al. (2001) observed the near- infrared CO-band 
head lines, where dust obscuration is far less of a concern 
than in the optical. Their observed sigma-drops, of gener- 
ally < 20 per cent, were due to nuclear bars or discs, vis- 
ible in part because of their sub-arcsecond seeing. Chung 
& Bureau (2004) subsequently reported on the detection of 
sigma-drops in a larger sample of edge-on galaxies. They too 
found decrements of up to 15-20 per cent due to likely cen- 
tral discs. These papers revealed that the sigma-drops are 
only present in a small fraction of galaxies, and that they 
provide a smaller perturbation compared to that capable of 
a large scale bar. 

4.3 A first order correction to a for bar dynamics 

As already discussed, we found that the 20 galaxies with 
large scale bars, from among the 64 predominantly inactive 
galaxies with direct SMBH mass measurements, are some- 
what offset from the barless Mbh-o" relation towards higher 
velocity dispersions. In this section we discuss how this may 
be linked to orbital bar dynamics. We do not, however, con- 
sider the potential bias due to nuclear/inner bars, or nuclear 
discs with young stars which can exist in both barred and 
non-barred galaxies. In what follows we offer a simplistic, 
statistical correction for the large scale bars. 

When measuring the velocity dispersion in the inner- 
most portion of an edge-on SA galaxy, one has contribu- 
tions mainly from the bulge orbits. This is because in the 
zero-th order approximation, the orbits of disc stars are cir- 
cular. Consequently, any disc stars seen (in projection) at 
the very center of the galaxy will have velocities perpendic- 
ular to our line-of-sight. Of course the disc star orbits are 
not exactly circular, and the measurements refer not only 
to the very center but to a region around it (see the previ- 
ous subsection). Inspite of this, the contribution of the disc 
to the velocity dispersion will be small compared to that of 
the bulge and can be neglected. This is not the case for SB 
galax;ies. 

The situation is more complex in a barred galaxy. Here 
the building blocks are closed periodic orbits, whose shape 
changes a lot as a function of location within or around the 
bar (e.g. Athanassoula 1992, her figure 3). If these are sta- 
ble, they trap around them regular orbits, which contribute 
to the velocity dispersion, in the same way as the orbits 
trapped around the near-circular orbits in SA galaxies. In 
the SB galaxies, however, there are further effects. Since the 
orientation of the bar is random with respect to the line of 
nodes, the velocity along the orbit is not necessarily perpen- 
dicular to the line-of-sight unless the bar is along the line 



of nodes, or is perpendicular to it (but, in this second case, 
for a much smaller inner radial range). Furthermore, some 
periodic orbits have loops which, if projected on the centre, 
can further increase the velocity dispersion. Readers wishing 
to understand these effects better can find more extensive 
explanations in Bureau & Athanassoula (1999), where indi- 
vidual periodic orbits are analysed. The bottom line is that 
the orbital structure of the bar will increase the observed 
velocity dispersion beyond that of the bulge alone, and will 
thus partly explain the (positive) sign of the deviation of 
barred galaxies from the regression line obtained from the 
remaining sample. Furthermore, the simulations analysed 
in Bureau & Athanassoula (2005) show that these effects 
can be quite strong and may enhance the observed velocity 
dispersion by as much as 10 - 40 per cent, depending on 
the bar's strength and orientation (Bureau & Athanassoula 
2005, their figure 1). This range is in agreement with the 
amplitude of the observed deviations. It is thus clear that 
these deviations may well be due to the orbital structure 
characteristics due to the bar potential. 

Here we provide a simplistic estimate of the affect from 
large-scale bars on the velocity dispersion (in the disc plane). 
Imagine a bar with major- and minor-axis lengths of 4 and 
2 kpc, respectively. The bar thus has b/a — 0.5 and is some 
8 kpc long. Seen end on, the bar is 4 kpc wide. To give some 
indication as to what fraction of the bar may be sampled by 
apertures used to measure a galaxy's central velocity disper- 
sion, 1" corresponds to ~0.1 kpc at a distance of 20 Mpc. 
Central velocity dispersions are typically measured within 
the inner couple of arcseconds, thus sampling the inner ±0.1 
kpc in this example. The main family of orbits sustaining the 
bar are xl orbits, which are like ellipses that are aligned with 
the bar and centred on the bar; these are of course different 
from purely radial orbits. As noted above, when this bar is 
viewed end-on, and at 90 degrees from this orientation (i.e. 
viewing the full length of the bar) , the central stellar orbits 
will again be perpendicular to our line of sight. This, how- 
ever, will not be the case at intermediate viewing angles. 
That is, the increase in velocity dispersion will be due to 
elliptical orbits entering the 0.1-0.2 kpc aperture when the 
bar has a position angle different from or 90 degrees, with 
the maximum when the bar is at 45 degrees, contary to what 
one may at first imagine. 

To roughly correct the observed central velocity disper- 
sion (cTobs) for broadening by non-circular disc motions due 
to a large-scale bar, and thus acquire the bulge velocity dis- 
persion (cTbuigc), may require the use of an expression like 



bulge *^obs 



CKotSin^(i) X sin {PA) x cos^{PA), (11) 



where i is the inclination of the disc such that i = repre- 
sents a face-on orientation, and C is a value between and 1 
and proportional to the bar's strength. The position angle of 
the bar, relative to the disc's apparent minor axis, is denoted 
by PA. For the reason discussed above, the above correction 
is minimal when the bar is parallel to the projection axes 
(either major or minor axis). In passing we note the compli- 
cation that as a disc's inclination increases from face-on (0 
degrees) to edge-on (90 degrees) , the observed position angle 
of the bar can appear increasingly aligned with the (appar- 
ent) position angle of the disc (e.g. Debattista et al. 2002). 
Measured disc inclinations themselves can also be subject 
to the thickness of the disc, departures from circularity (e.g. 
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Table 3. Disc rotations and inclinations, and bar orientations 



Gal. Id. 



Type 



50 



100 

0" 



bulge 



200 

[km s ] 



Figure 7. Adaptation of Figure [5] such that (i) the central veloc- 
ity dispersions of the barred galaxies, cr ohs^ have been adjusted to 
a new "o"^J^lgJ," value following eguation ll ll using C = 0.7, and (ii) 
the lines and shaded area now trace the A/bh-f relation defined 
by the non-barred galaxies, as given by the second last entry in 
Table [2] Squares, stars and circles denote elliptical, barred and 
non-barred disc galaxies, respectively. The small lines emanating 
from the barred galaxies shows their previous location in Figure[5] 



Andersen et al. 2001) and potential ellipticity gradients in 
which the outer isophote shape depends on exposure depth. 

In an effort to test if bar dynamics may be a viable 
(partial) explanation for the offset barred galaxies in the 
Mbh-o" diagram, we have, for an initial rough investigation, 
taken the disc inclination, position angle and Kot measure- 
ments from Hyper leda (see Table [3]), and assumed a C- value 
of 0.7 for all. We note that the actual C- value is currently 
unknown, and likely to be different for each galaxy, and we 
caution that the heterogeneous nature of the data that goes 
into HyperLeda will also be a source of scatter. Given the 
simplistic nature of equation 1111 we were content, as a first 
approach, to estimate the bar's position angle by eye. Fig- 
ure [7| shows how the barred galaxies shift when using our 
rough estimate of (Jbuigc rather than Uobs. One can see that, 
as expected, the corrections are significant only in the cases 
where Uobs, and thus }Abh, are small, and that it shifts these 
cases towards smaller (Tbuigc values. 

Depending on the values of C, this correction has the 
potential to improve the correlation. Obviously though, this 
work requires greater investigation than we are able to afford 
here. We also again remark that we have not dealt with 
the second, independent phenomena known as sigma-drops, 
which may potentially even account for some of the apparent 
over-shoot using equation [TT] 



4.4 Nuclear star clusters 

The offset nature, or rather the location below the upper 
envelope of points, of some of the apparently non-barred 
AGN in the A/bh-o" diagram (Figure [S]:) prompts one to 
consider alternative scenarios to bar dynamics. A couple of 
non-barred galaxies with direct SMBH mass measurements 
also appear to be offset from the barless Mbh-o" relation (see 
Figure [2]). If they are not pseudo bulges, in which the bar 
has disappeared, then something other than elevated veloc- 
ity dispersions may be responsible for their deviant nature. 
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The maximum disc rotation velocities, Vmax, have been taken 
from HyperLeda. The disc inclinations have been obtained from 
available photometry, under the assumption that any intrinsic 
disc ellipticity is small and that the observed minor-to-major axis 
ratio is due to disc inclination (* taken from Graham &; Li 2009). 
The rough position angle of the bars, P.A.ban relative to the 
minor axis, has been estimated by eye. 



It has already been noted by Graham (2008a) that the 
inactive galaxy NGC 2778 has only a weak bar, yet this 
nucleated galaxy is below the barless Mbh-o" relation by 0.7 
dex. Curiously, the mass of the nuclear cluster in this galaxy 
is '^S times (i.e. 0.7 dex) greater than the mass of its black 
hole. Furthermore, the Milky Way and M32 have nuclear 
star cluster masses, Mnc, which are lOx more massive than 
their SMBH masses (Graham & Spitler 2009, and references 
therein) . 

Nuclear star clusters (Ferrarese et al. 2006; Wehner & 
Harris 2006) are now known to coexist with SMBHs at the 
centres of low- and intermediate-mass spheroids (Filippenko 
& Ho 2003; Graham & Driver 2007a; Gonzalez Delgado et 
al. 2008; Seth et al. 2008; Graham & Spitler 2009). The 
star clusters undoubtedly contribute, at some level, to the 
feeding of the central black hole, whether by direct stellar 
infall (e.g., Lightman & Shapiro 1978; Merritt & Vasiliev 
2010) or through stellar winds (e.g., Ciotti et al. 1991; Soria 
et al. 2006; Hueyotl-Zahuantitla et al. 2010). Schartmann et 
al. (2010), for example, describe how stellar mass loss from 
the nuclear star cluster (or nuclear disc, Davies et al. 2007) 
in NGC 1068 can feed the black hole of this active Seyfert 
galaxy. In unrelated work, Bekki & Graham (2010) argue 
why a seed SMBH in nuclear star clusters may be a necessary 
ingredient to explain the absence — through SMBH binary 
heating and ultimately erosion — of nuclear star clusters in 
massive galaxies built via hierarchical merging. Given the 
growing number of probable connections between SMBHs 
and nuclear star clusters, which is no doubt yet to be fully 
appreciated, it seems pertinent to explore the inclusion of 
the nuclear star cluster. 

While we do not have nuclear star cluster masses for the 
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AGN sample with (reverberation-mapping)-derived SMBH 
masses, and can not therefore re-derive the /-factor, progress 
has commenced on acquiring star cluster masses for the lo- 
cal sample of predominantly inactive galaxies with direct 
SMBH mass measurements (Graham & Spitler 2009). In ad- 
dition, Graham (2010, in prep.) tabulates yet more galaxies 
having both a massive black hole and a nuclear star clus- 
ter. Following Graham & Guzman (2003) and Balcells et al. 
(2003, 2007), the nuclear star cluster fluxes have been de- 
rived by simultaneulsy fitting for the cluster, the host bulge 
and, when present, nuclear and large scale discs. Failing to 
account for these components can bias the Sersic model used 
to describe the bulge light, and thus bias the flux of the 
nuclear star cluster. Foreground stars, background galaxies, 
and any apparent internal dust is masked out before the see- 
ing convolved models are fit to each galaxy's distribution of 
stellar light. The nuclear cluster fluxes are then converted 
into masses following the methodology in Graham & Spitler 
(2009). 

With this expanded data set, Figure[H]explores the spec- 
ulation by Graham & Spitler (2009) that an insightful quan- 
tity to plot on the vertical axis may be the combination of 
the SMBH plus the nuclear star cluster mass. Nayakshin et 
al. (2009) discuss how competing feedback from a SMBH 
and a nuclear star cluster may explain either end of such 
an (Mbh + Mnc)-<y relation (see also McLaughlin, King & 
Nayakshin 2006 and Cantalupo 2010). When both types of 
nuclei exist, it may make sense to combine their masses, and 
Figure [S] presents the first ever "{M+M)-a" diagram show- 
ing this. Rather than treating galaxies as if they only have 
one type of nuclear component. Figure |8] shows that the low 
mass end of the Mbh-cr diagram flattens when transformed 
into the (Mbh + -A/nc)-o" diagram. While too far off topic 
for the present paper, Graham (2010, in prep.) expands this 
analyis to include (the transition to) low-mass galaxies for 
which only the nuclear star cluster mass is known. 

If NCs have a different velocity dispersion than their 
host bulge, as is the case for NGC 205 (Carter & Sadler 
1990, see their figure 2), then they could bias measure- 
ments of the host bulge's central velocity dispersion a. How- 
ever, for most galaxies the NC-to-(host bulge) flux ratio will 
be small in typical ground-based aperture measurements, 
which could explain why observed "sigma-drops" are rela- 
tively rare, which is not to say that they do actually exist in 
great numbers. High spatial resolution spectra would be de- 
sirable for pursuing this issue of contamination at distances 
of the Virgo and Fornax galaxy clusters. 

It is interesting that the galaxies with the smallest 
SMBHs in the left hand side of Figure [8] tend to be dusty 
Seyfert 2 galaxies with AGN flux that contaminates the 
nuclear stellar fiux, and thus any NC mass measurement. 
These AGN with direct SMBH mass measurements are, like 
the reverberation-mapped AGN, ofi^set below (or rightward 
of) the upper envelope of points in the Mbh-o" diagram. 
Their AGN flux, and dusty central region, may hide a nu- 
clear star cluster and effectively bias these points low in 
the (Mbh + Mnc)-o" diagram. Nonetheless, from Figure [5] 
it is apparent that the slope at the low-mass end of the 
(Mbh+Mnc)-cr diagram is reduced relative to the Mbh-o" rela- 
tion. As discussed in Graham (2010, in prep.), this behavior 
appears to track the change in slope seen in the luminosity- 
(velocity dispersion) diagram (e.g. De Rijcke et al. 2005; 
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Figure 8. Adaptation of Figure [5] using the 64 galaxies with di- 
rect SMBH masses (plus NGC 4395, lower left). The solid upward 
pointing arrow heads show the location of the combined black hole 
masses (Graham 2008b and Table [1} plus the nuclear star clus- 
ter mass (Graham & Spitler 2009 and Graham 2010, in prep.). 
For seven systems with only an upper limit to the nuclear cluster 
mass, a bi-directional arrow denotes the possible (Mbh + Mnc) 
mass range. Galaxies with dust obscured nuclei, and/or signifi- 
cant AGN whose flux may hide a nuclear star cluster, are plotted 
with open triangles on top of the mass of the SMBH. 



Matkovic & Guzman 2005, and references therein), and is 
thus connected with the (A/bh + j^fnc)/-f'buigc ratio. 



4.5 Radiation Pressure 

Marconi et al. (2008) explored the effect that radiation pres- 
sure (e.g. Mathews 1993), from the AGN's continuum emis- 
sion, may be having on the surrounding gas clouds, and thus 
the effect on measurements of the virial product. Such an 
outward force would counteract, at some level, the inward 
force of gravity on the gas. Marconi et al. (2008) therefore ar- 
gued that neglecting the radiation pressure in reverberation 
mapping analyses results in a systematic under-estimation 
of true black hole masses. Armed with flux measurements 
at 5100A, they empirically derived an /-factor of 3.lli'5, 
smaller than the value of 5.5±1.7 in Onken et al. (2004), 
but now with an extra luminosity-dependent term. 

Marconi et al.'s (2008) radiation pressure term will re- 
sult in an upward revision to the masses of high-accretion 
AGNs. Given that Narrow Line Seyfert 1 (NLSl) galaxies 
tend to have high accretion rates, and may reside below the 
Mbh - o- relation (Wandel 2002; Mathur & Grupe 2005; Wu 
2009; but see Botte et al. 2005; Komossa & Xu 2007; Decarli 
et al. 2008a), Marconi et al. (2008) point out that the ap- 
parent offset of the NLSl galaxies relative to the broad- line 
AGN may be due to past under-estimates of their black hole 
masses. 

While Netzer (2009) concluded that radiation pressure 
is not significant, and that the empirically-determined value 
by Marconi et al. (2008) is too high, the rebuttal by Marconi 
et al. (2009) counters this. What remains unclear is whether 
the slight reduction in the scatter for the AGN Mbh — "" 
relation is due to an improvement in the physical model 
or simply the addition of an extra fitting parameter (see 
Peterson 2010). Here we simply remark that a corrective 
term for radiation pressure would act to drive the /-factor 
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to even lower values. That is, while our present work has 
halved the commonly used value, it may still need to be 
revised further downward. 



4.6 Implications for AGN 

As a measure of the relationship between emission line 
widths and the gravitational potential of the central SMBH, 
the /-factor is a direct constraint on the kinematics and ge- 
ometry of the broad-line region. An / value of 3 is implied 
for the simplest model of a BLR (Netzer 1990). The next 
level of sophistication in BLR modeling involves a flatten- 
ing of the spatial distribution. Thin discs, thickened discs, 
or composites of planar and random motions have all been 
modeled by various groups (see discussions by Krolik 2001; 
CoUin et al. 2006; Labita et al. 2006; Decarh et al. 2008b). 
The uniqueness of any of these solutions remains unclear, 
but improved estimation of both / and the scatter around / 
(to be expected with an inclination-dependent BLR model) 
are positive steps towards understanding the BLR geometry. 

One consequence of our downward revision in the mean 
value of / is to reduce the black hole masses for any AGNs 
calibrated to the Onken et al. (2004) result. Perhaps frus- 
tratingly, our lower value of / gives masses quite similar 
to those which would have been produced by the simplis- 
tic BLR model which had been in use prior to the work by 
Onken et al. (modulo the relationship between FWHM and 
ciino that one adopts). Obviously, the change in / directly 
impacts any reverberation results that adopted the Onken 
et al. (2004) value (Peterson et al. 2004, 2005; Metzroth et 
al. 2006; Bentz et al. 2006, 2007, 2009a; Denney et al. 2006, 
2009; Sergeev et al. 2007; Grier et al. 2008). 

Furthermore, the single-epoch methods of SMBH mass 
estimation which were calibrated to the reverberation map- 
ping results would be similarly affected (Vestergaard & Pe- 
terson 2006; McGiU et al. 2008; Wang et al. 2009; Greene, 
Peng & Ludwig 2010). Because these single-epoch methods 
are the primary basis for estimating black hole masses from 
AGN spectra, the effects of changing / can be far-ranging, 
from individual objects to large AGN surveys (Barth et al. 
2005; Woo et al. 2006, 2008; Dong et al. 2007; Zhang et al. 
2007; Treu et al. 2007; Vestergaard et al. 2008; Shen et al. 
2008a,b; Trump et al. 2009; Merloni et al. 2010; Greene et 
al. 2010; Morgan et al. 2010; Lamastra et al. 2010). 

For some applications, a new / value has no significant 
ramifications: e.g., the evolution in the relationship between 
black holes and spheroids in AGN samples (e.g. Peng et al. 
2006a,b; Bennert et al. 2010; Decarh et al. 2010a,b); or the 
slope of the active BH mass function (e.g. Kelly, Vestergaard 
& Fan 2009; Vestergaard & Osmer 2009). In other cases, 
the fact that the modifications at present are uniform to 
all BH masses means that the results are simply shifted by 
a factor of ~ 2 (e.g., the Eddington ratio distributions of 
AGNs: KoUmeier et al. 2006; the position of the peak in 
the black hole mass function: Kelly et al. 2009). In some 
circumstances, the shift in black hole masses may have a 
still larger impact (e.g., in modeling the growth of black 
holes over cosmic time: Yu & Lu 2008; Vestergaard & Osmer 
2009). 

However, it bears noting that, if / differs for barred and 
non-barred AGNs, changes in the bar fraction with redshift 
and/or with galaxy mass could dramatically alter the con- 



clusions of studies unaffected by a simple / shift (i.e., the 
slope of the BH mass function, and the evolution of Mbh- 
galaxy relationships). 

As a final example of potential implications, we dis- 
cuss the growing evidence that massive black holes in high- 
redshift AGNs may pose a challenge to the notion that 
such objects can grow from stellar-mass seeds via normal 
(Eddington-limited) accretion in the time available since the 
Big Bang (Dietrich & Hamann 2004). To explore the impact 
of a revised / value on such analyses, we focus on two par- 
ticular objects. First, the 3 x 10® Mq black hole found in 
the z = 6.41 AGN, SDSS J1148-h5251 (Willott, McLure & 
Jarvis 2003), was measured by a mass equation which as- 
sumed f = 1 (for a line width measured by the FWHM). 
Our / value would increase the BH mass for this object, 
making it even harder to grow such a massive BH at such 
an early time. In contrast, the BH mass in the z — 2.131 
AGN, Q0019H-0107, of 9.5x10® Mq (Dietrich et al. 2009) 
would be halved, and would therefore no longer require a 
seed black hole quite so massive as 10"" Mq. These two 
cases demonstrate the complexity involved in comparing any 
SMBH masses in AGNs from the literature, especially since 
the pedigrees of such estimates are sometimes not clearly 
documented. 



5 CONCLUSIONS 

Using a sample of 64 galaxies with directly measured su- 
permassive black hole masses, the "classical" or "standard" 
Mhh-<y relation (which contains galaxies of all morphological 
type) is shown here to have a steeper slope than previously 
recognised. Due to (i) the inclusion of more barred galaxies, 
which tend to be offset from the non-barred galaxies (Gra- 
ham 2008a,b) and had previously been under-represented 
in the Mbh-f diagram, along with (ii) increased black hole 
masses at the high-mass end of the Afbh-o" diagram (e.g. 
Gebhardt & Thomas 2009; Shen & Gebhardt 2009; van 
den Bosch & de Zeeuw 2010), and (iii) the use of a 10 
per cent, rather than a rather optimistic 5 per cent, un- 
certainty on the velocity dispersion, the slope is now some- 
where between 5 and 6 depending on the type of regres- 
sion used (see Table [2]). From a non-symmetrical regres- 
sion of logA^bh on log a, and using a 10 per cent uncer- 
tainty on the velocity dispersion, one has log(Afbh/AfQ) — 
(8.13±0.05)-h(5.13±0.34)log[t7/200kms"^], withanr.m.s. 
scatter of 0.43 dex. While, formally, the intrinsic scatter is 
0.32 dex, we note that a more detailed analaysis of the un- 
certainty in the velocity dispersions would be welcome. 

If intermediate mass black holes exist (e.g. Greene & Ho 
2004, 2007), and an observational selection bias is artificially 
truncating data in the Afbh-o" diagram below Mbh = 10® Mq, 
then the above relation is biased. Using a regression of 
log a on log Afbh gives the bias- free relation log(Afbh/A/Q) = 
(8.15 ± 0.06) + (5.95 ± 0.44) log[CT/200kms"^]. When one 
doesn't know the morphological type of their galaxy, these 
expressions are currently the best available estimators for 
Mhh from measurements of a. 

As was first pointed out by Graham (2007a, 2008a, b) 
and Hu (2008), the elliptical-only galaxies, and the non- 
barred galaxies, define tighter relations with less scatter and 
a reduced slope than is obtained when using the full galaxy 
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sample. Furthermore, the barred Mbh-o" relation is shown 
here to reside 0.47 and 0.45 dex below the elliptical-only 
and barless Mbh-o" relations, respectively. When one knows 
the morphological type of their galaxy, these relations (given 
in Table [2]) are preferred for predicting Mbh from measure- 
ments of a (as described in Section [2.1. 2 1) . The barless and 
elliptical-only Mbh-o" relations have a total r.m.s. scatter of 
0.37 and 0.34 dex, respectively, and a slope of 4.95 and 4.87 
when constructed with a symmetrical treatment of the data. 
A slope which is also consistent (at the 1-sigma level) with 
a value of 5 is obtained from a linear regression of log a on 
log Mbh (see Tabled]). 

In terms of an integer slope for the Afbh-o" relation, when 
constructed to address the "Theorists Question" (Novak et 
al. 2006), a value of 5 is thus preferred for elliptical galaxies, 
in agreement with the prediction by Silk & Rees (1998). 
Due to both new and refined SMBH masses, and increased 
values at the high-mass end, this slope has increased from 
the values around 4 that were reported by Graham (2008b), 
Hu (2008) and Giiltekin et al. (2009b). The even steeper 
slope, close to a value of 6, which is obtained for the full 
galaxy sample is due to the relatively offset nature of the 
barred galaxies in the lower-left of the Mbh-o" diagram and 
the use of a linear regression not biased by sample selection 
which excludes BHs with Mbh < W'^Mq. 

Using a sample of 28 AGN with available virial products 
and host bulge velocity dispersion measurements, we have 
explored their location in the A/bh-o" diagram. We have then 
derived a new /-factor for converting AGN virial products 
(equation [5]) into black hole masses. Our value / = 2.8lg'5 
is a factor of two smaller than the commonly used value 
5.5 ± 1.7 (Onken et al. 2004, see also Woo et al. 2010 who 
reported / — 5.2 ± 1.2). Moreover, this value might come 
down even further due to processes that may be acting to 
enhance the offset nature of AGN virial products from the 
Mbh-o" relation defined by predominantly inactive galaxies 
with direct measurements of their supermassive black hole. 
For example, radiation pressure from the AGN may be act- 
ing to partly counter-balance the inward force of gravity (on 
the broad line region clouds) due to the black hole, possibly 
resulting in anemic virial products (Marconi et al. 2008). 
Correcting for this possibility will only lower the derived 
/-factor further. As such we conclude that our reduced /- 
factor of / = 2.8'ta:l (= 2.6+^j when excluding NGC 5548) 
may still be an upper limit, and therefore some quasar and 
AGN masses in the literature are too high by at least a fac- 
tor of 2. Some of the implications of this are discussed in 
Section 1321 

This paper has also highlighted a number of issues. 
First, the impact of galaxy type (i.e. elliptical, barred, non- 
barred) and thus sample selection on the A'lhh-cr relation is 
important. Depending on the relative numbers of barred and 
non-barred galaxies, the best-fitting Afbh-o" relation changes. 
This of course undermines past efforts which had focussed 
on the exact value of the slope. It also probably voids several 
evolutionary studies based on earlier Mbh-o" relations. The 
present study is also not immune from this, although the 
construction of barred and non-barred Mbh-o" relations for 
galaxies with directly measured SMBH masses is a positive 
step, plus this study is based on a larger and more repre- 
sentative sample of galaxies than ever before. Second, as the 
black hole masses continue to be updated and refined, sys- 



tematic shifts in their masses can affect the Mbh-o" relation. 
For example, van den Bosch & de Zeeuw (2010) have shown 
how past assumptions that galaxies are either axisymmetric 
oblate spheroids or spherical, rather triaxial, has resulted 
in an underestimation of their SMBH mass. Gebhardt & 
Thomas (2009) have revealed how ignoring the dark mat- 
ter's infiucnce on the stellar dynamics can bias the stellar 
mass-to-light ratio resulting in an underestimation of the 
SMBH mass. In the unlikely event that all of the current 
black hole masses at the high mass end are doubled, i.e. in- 
creased by 0.3 dex, it would result in a 10 per cent increase 
to the slope of the relation from 10^ to 10^ Mq (i.e. over a 
range of 3 dex in SMBH mass). Third, this work confirms 
Graham's (2007a;2008a,b) and Hu's (2008) finding that the 
scatter about the Afbh-o" relation is quite large, and that ve- 
locity dispersion alone is therefore probably not the driving 
force which dictates SMBH mass; if it was all that mattered, 
then the barred galaxies would follow the same distribution 
as the non-barred galaxies in the Mbh-o" diagram. 

There are several potential causes for this factor of ~3 
offset between the barred and non-barred galaxies in the 
Mbh-o" diagram, although it remains unclear what their var- 
ious contributions are. Possibilities include heightened ve- 
locity dispersion measurements due to (i) elongated motions 
along a bar, (ii) vertical instabilities taking stars out of the 
disc plane, (iii) rotational shear due to approaching and re- 
ceding parts of a disc within one's aperture, and iv) the 
variety of orbital shapes in and around the bar. The previ- 
ously ignored contribution from massive nuclear star clusters 
may also play a role (see Section [4. 2[) . Nuclear star clusters 
are also prevalent in low-mass, elliptical galaxies. Including 
more of this galaxy type, to check for their offset behaviour 
in the Mbh-o" diagram, may be insightful. As a first step, we 
have presented the first ever (Mbh-l-Mnc)-o" diagram (Fig- 
ure [8|, qualitatively showing how the slope at the low-mass 
end of the {M+M)-a diagram is shallower than it is at the 
high-mass end where nuclear star clusters are not detected. 
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After completion of this work we became aware of four 
additional velocity dispersion measurements for our AGN 
sample. Using a direct filtering method, Garcia-Rissmann 
et al. (2005) report velocities of 116±20, 153±24, 83±11 
and 125±12 km s"^ for NGC 3783, 4593, 6814 and 7469, 
respectively. Within the quoted errors, these values agree 
with those used in this paper. 
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Table 4. AGN data 



Gal. Id. 


Class. 


Virial Product 


fj 






[lO^Mo] 


[km s-i] 


1 


2 


3 


4 






28 non-barred 


galaxies 


3C 120 


SO (pec?) 


lo.iti;! 


162+20^ 


Ark 120 


Sb pec 


27.2+3.5 


239+36", 221 ± 17^, 224 ± 15 


Arp 151 


SO pec 


1 99+0. 17a 


124+120,118 + 4^, 119 + 4 


Mrk 110 


Sc 


4.56+1.08 


86±13'3,91±25^, 87 ±12 


Mrk 202 


Sc 


0-26+g:n^ 


86+14^,78 + 3^, 78 ±3 


Mrk 279 


SO 


6.35+1.66 


197+12^ 


Mrk 290 


El 


4 42+0.67B 
^■^^-0.67 




Mrk 335 


so/a 


2.58+0.67 




Mrk 509 


E2 


26.1+2.1 




Mrk 590 


SA(s)a 


8.64+1.34 


194+20*^, 189 ± 6^, 169 ± 28^, 189 ± 6 


Mrk 1310 


E 


r|4i+0.16A 


50+16^,84 + 5^, 81 ±5 


NGC 4395 


Sm 


0.065±0.020<^ 


25-35-f*, 30 ±2.5 


NGC 4748 


Sa 


Q47+o.igA 

^■^' -0.23 


105+13^ 


NGC 5548 


(R)SA(s)0/a 


-,9 -,+0.5 J 


183±27<3, 201 ± 12^^, 195 ± 13^, 197 ± 8 


PG 0026+129 


El 


71.4+17.4 




PG 0052+251 


Sb 


67.3+13.7 




PG 0804+761 


E3 


126.4+14.5 




PG 0844+349 


Sa 


16.8+7.0 




PG 0953+414 


E4 


50.1+10.7 




PG 1226+023 


E3 


160.6+34.1 




PG 1307+085 


E2 


80.1+22.1 




PG 1411+442 


E4 


80.5+26.5 




PG 1426+015 


E2 


236.7+69.8 


217+15^ 


PG 1613+658 


E2 


50.7tij 




PG 1617+175 


E2 


108.0+25.3 


183+47^ 


PG 1700+518 


El 


-,49+33 




PG 2130+099 


(R)Sa 


6.9+2.7° 


172+46^ 


SBS 1116+583A 


Sc 


1.05t°:g^ 


50+18^,92 + 4^, 90 ±4 






14 barred galaxies 


Fairall 9 


SBa 


46.3+10.0 


228+20^ 


Mrk 79 


SBb 


9.52+2.61 


130±20'3, 130 ± 12^^, 130 ± 10 


Mrk 142 


SBO/a 


0.40«:1*^ 




Mrk 817 


SBc 


9.46+1.24^ 


142±21<3, 120 ± 15^^, 127 ± 12 


NGC 3227 


SAB(s) pec 


1 00+0.29 s 

j^"5y_o.3i 


131+11", 136 ± 4^^, 128 ± 13^, 135 ± 4 


NGC 3516 


(R)SB(s) 


O-'O-0.76 


144+35^, 164 ± 35^, 181 ± 5^, 180 ± 5 


NGC 3783 


(R)SB(r)a 


5.41+0.99 


95+10" 


NGC 4051 


SAB(rs)bc 


r, O1+0.10B 

"^■■'^-o.og 


84+9*5, 89 ±3^, 88 ±13^, 88 ± 3 


NGC 4151 


(R)SAB(rs)ab 


o O-I+1.04F 


93+14*3,97 + 3^,119 + 26^, 97 ± 3 


NGC 4253 


SBa 


0.32«-J«^ 


93+32^ 


NGC 4593 


(R)SB(rs)b 


i.8j:«:l« 


124+29^, 135 ± 6^, 135 ± 6 


NGC 6814 


SBbe 


o •5fi+0.63A 
•3-JO_o.64 


115+18^,95 + 3^, 96 ±3 


NGC 7469 


(R)SAB(rs)a 


2.21+0.25 


152+16", 131 ± 5^^, 133 ± 5 


PG 1229+204 


SBc 


13.4+6.2 


162+32^ 






4 excluded gi 


ilaxies 


IE 0754.6+3928 


7 


91 n+40.3H 




3C 390.3 


Sa 


52.3+11.7 


240+36^', 273 ±16^, 268 ± 15 


IC 4329A 


SAO 


1.80_2 16 


122+13" 


PG 1211+143 


E2 


26.5+7.9^'^ 





Column 1. Galaxy Identification. Column 2. Galaxy Classification from NED, except for Fairall 9 and PG 1229+204 which have been 

denoted as barred by Bentz et al. (2009b, their Table 5), plus Mrk 202 and SBS 1116+583A whose type has been assigned based on 

SDSS images and colour. Column 3. Virial Product (and also black hole mass if / = 1) using the data from Peterson et al. (2004) 

unless otherwise noted. Column 4. Stellar velocity dispersions; the bold value is the weighted mean. 
References and notes: ^ Bentz et al. (2009a), HP line. ^ Denney et al. (2010). ^ Peterson et al. (2005). ° Crier et al. (2008). ^ 
Weighted mean from Peterson et al. (2004) and Denney et al. (2010). ^ Bentz et al. (2006). ^ Denney et al. (2006). ^ Sergeev et al. 
(2007). ^ Time-lag uncertain (Peterson et al. 2004) and thus Virial Product not used. '' Weighted mean from Peterson et al. (2004), 
Bentz et al. (2009a) and Denney et al. (2010). ^ H^ 
(2004). ^ Nelson et al. (2004). ° Greene & Ho (2006). 
uncertainty. ^ Filippenko & Ho (2003, their section 3). 



and H^, but not Hq,, line used here. ^ Nelson & Whittle (1995). ^^ Onken et al. 



et al. 



^ Woo et al. (2010). '^ Ferrarese et al. (2001), using their suggested 15 per cent 
S Watson et al. (2008). ^ Dasyra et al. (2007). ^ Oliva et al. (1995). ^ Arribas 
(1997). ^ Green et al. (2004). 
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